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Abstract 

Previous work in the literature had built a formalism for spatially averaged equations for the 
scale factor, giving rise to an averaged Raychaudhuri equation and averaged Hamiltonian con- 
straint, which involve a backreaction source term. The present paper extends these equations to 
include models with variable Newton parameter and variable cosmological term, motivated by the 
nonperturbative renormalization program for quantum gravity based upon the Einstein - Hilbert ac- 
tion. We focus on the Brans -Dicke form of the renormalization - group improved action functional. 
The coupling between backreaction and spatially averaged three-dimensional scalar curvature is 
found to survive, and a variable-G cosmic quintet is found to emerge. Interestingly, under suitable 
assumptions, an approximate solution can be found where the early universe tends to a FLRW 
model, while keeping track of the original inhomogeneities through three effective fluids. The re- 
sulting qualitative picture is that of a universe consisting of baryons only, while inhomogeneities 
average out to give rise to the full dark-side phenomenology. 
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I. INTRODUCTION 

Over the last ten years, the use of the effective average action [l| and of the renormal- 
ization - group equations has made it possible to obtain encouraging evidence in favour of 



with all run- 



Einstein's gravity being renormalizable at non - perturbative level [2|, 
ning couplings having a finite limit at large momenta. The cornerstone of this program is 
the discovery of a new non- Gaussian ultraviolet fixed point, besides the trivial one at the 

origin aaaa 

. The resulting picture seems to be as follows : at sub - Planckian distances, 
spacetime is a fractal. It can be thought of as a self- similar hierarchy of superimposed 
Riemannian manifolds of any curvature. As one considers larger length scales where the 
renormalization- group running of the gravitational parameters comes to a halt, the space- 
time ripples gradually disappear, and a classical four - dimensional spacetime manifold is 
recovered Q]. 

In the simplest possible terms, the renormalization group improvement consists in the 
modified Einstein equations [l^ 

R^.. - Ig^^uR + Hk)9>^u = 87rG(A;)T^„ (1.1) 

where the Newton parameter G and cosmological term A are now dependent on the scale 
k, k being the running cut-off of the renormalization group equation |1|. In cosmology, the 
dynamical evolution is determined by a set of renormalization group equations by means 
of the cut-off identification k = k{t) which relates the energ y s cale of the running cutoff 



k of the renormalization group, with the cosmic time t. In [ll| it has been shown that, 
in a cosmological setting, the correct cutoff identification is oc t^^; it is thus possible 
to determine G{k{t)) and A{k{t)) in Eq. (1.1) once a renormalization - group trajectory is 
determined. 

At a deeper level, one integrates out all fluctuations with momenta larger than a cutoff k, 
and one takes them into account by means of a modified dynamics for remaining fiuctuation 
modes with momenta smaller than k. This modified dynamics is generated by a scale - 
dependent effective action r^, whose /c - dependence is ruled by the exact renormalization - 
group equations. Flow equations can be used for a complete quantization of fundamental 
theories. On denoting by S their classical action, one imposes the initial condition r{k = 
k) = S a.t the scale of the ultraviolet cutoff k, and one exploits the renormalization -group 
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equation to evaluate this averaged effective action Tk for all k < n, and one then takes the 
limits k —>■ and /? ^ oo. In the case of fundamental theories, the continuum limit n ^ oo 
exists after having renormalized finitely many parameters in the action, and is evaluated 
at a non- Gaussian fixed point of the renormalization- group flow. Such a new fixed point 
replaces the Gaussian fixed point which, at least implicitly, underlies the construction of 
theories which are instead perturbatively renormalizable [121 ] . 



Notwithstanding its merits (see, however, the criticism expressed in 13j), the renormal- 
ization group approach has only been applied, so far, to a strictly homogenous and isotropic 
universe. As is well known, however, the matter distribution in the observable universe 
may be considered homogenously distributed only on large scales so that it is interesting 
to investigate how this can affect the renormalization -group equations. Needless to say, 
lacking a detailed description of the matter distribution, we can only consider a kind of 
average universe obtained by averaging out the inhomogeneities. Fortunately, over the last 
decade, progress has been made on the longstanc iing problem of how to average a general 
inhomogeneous model. In particular, the work in [145] considers an irrotational fluid motion 
with the associated Einstein equations (in c = 1 units) 

R/iu - ^gtiuR + ^giiu = ?>-nGpu^Uy. (1.2) 

A flow - orthogonal coordinate system = {t,X^) is chosen, i.e. Gaussian normal coordi- 
nates comoving with the fluid. On writing = /^(X^,t), one has = = (1,0,0,0) 
and Ufj, = ^ = (—1, 0, 0, 0), where t is proper time. The spacelike hypersurfaces of constant 
t are assumed to foliate the spacetime manifold, and their flrst fundamental form, i.e. the 
spatial 3 -metric gij, is used to deflne 



J{t,X') = ^detg,,. (1.3) 

The spatial averaging of a scalar fleld ip as a function of Lagrangian coordinates and time 
on an arbitrary compact portion D of the fluid is deflned by the volume integral (cf. [isl ) 

(^(t, X'))n I Ht, X')J d'X, (1.4) 

Yd jd 

where the volume Vd is obtained as 

Voit) = JSX. (1.5) 
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A key role in the formalism is played by a dimensionless effective scale factor 

in terms of which the averaged expansion rate takes the form (the dot being used for the 
total derivative with respect to t) 



(9>« = ^ = 3^^. (1.7) 

Yd ao 



In [ij], the scale factor obeys the averaged Raychaudhuri equation 

3— + ^nG{p)n - A = Qz. = ^((^ - {9)nf)D - 2{a')n, (1-8) 
an o 

where is the backreaction source term and o"^ is the rate of shear, and the Hamiltonian 
constraint 

3 _ SnG{p)o + le^R)D - A = -^Qn, (1.9) 

where '•^^i? is the scalar curvature of the constant time spacelike hypersurfaces used in the 
spacetime foliation. 

Our aim here is to generalize the above averaging procedure to the variable - G cosmolo- 
gies resulting from the renormalization - group approach, to gain a better understanding of 
the conditions under which a FLRW universe can be recovered from a renormalization - 
group approach (rather than imposing a FLRW symmetry as has been done so far). For 
this purpose, section 2 derives the Buchert average of the field equations obtained from a 
Brans -Dicke approach to the renormalization -group improved gravitational action. The 
integrability condition for the above equations is obtained in section 3, while a variable -G 
cosmic quintet is found to emerge in section 4. Section 5 studies accelerating patches and 
stationary models, while a solution formula for the spatially averaged scalar curvature is 
obtained in section 6. A particular solution of the averaged equations is then investigated 
in Section 7, where we also consider the case of a nearly homogenous universe through the 
useful introduction of a set of effective fluids clarifying the role of the different terms. Results 
and open problems are discussed in section 8. The basic identities used for the evaluation 
of Buchert averages are described in the appendix. 
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II. THE BUCHERT AVERAGING METHOD WITH VARIABLE G AND A IN A 
BRANS - DICKE APPROACH 



In light of equations (1.3)-(1.9), it is clear that the Buchert averaging method aims at 
taking spatial averages of equations obtained from suitable contractions and traces of tensor 
field equations. In particular, within the framework of running- G models, an interesting 
class is given by the Brans -Dicke approach developed in 16|. Here, the Einstein - Hilbert 
action is renormalization- group improved by replacing the Newton constant and the cosmo- 
logical constant by scalar functions G and A in the corresponding Lagrangian density 16|. 
The position dependence of G and A is governed by a renormalization -grou p e quation, and 
they have the status of externally prescribed background fields (whereas in 17| G obeys an 
Euler - Lagrange equation), while the metric satisfies an effective Einstein equation similar 
to that of Brans -Dicke theory l™. 

n 

Following [16|, we assume that the total action functional S of our theory consists of 
the Einstein - Hilbert part Seh plus matter Sm plus a term describing the four - momentum 
carried by the fields G{x) and A(a;), i.e. (hereafter (p denotes the collection of matter fields 
coupled to gravity) 

S = Seh{9, G, A) + Sm{9, + Seig, G, A). (2.1) 
For the explicit form of the terms in (2.1) we refer the reader to Ref. ^], for length reasons. 
The resulting renormalization -group improved Einstein equation is given by 



R 



fJ,U 



23t^uR 



SttGT, 



flU 



where t^^ results from the x dependence of G [16| : 

^ {2(V^G')(V.G) - GV^V.G - g,,[2{V ,G){VfG) - Gg^^-'V.V.G] } 



r, 



G^ 



(2.2) 



(2.3) 



while "dfj^i, is obtained from the functional derivative of Sg with respect to the metric, and 
can be taken to be of the form jl6| 



1 



iy^G)iy,G)--g,,iy,G)iy''G) 



(2.4) 



3 1 

'^^ ~ ~2G2 

which is a Brans -Dicke energy - momentum tensor for the field 1/G. For our purposes it is 
convenient to add explicitly the tensors r^i, and 'dfj,u to find 

= ^ {^(V^G)(V.G) - GV^V.G + g,, [-^(VpG)(V''G) + G\3g\ } , (2.5) 



where Q is the standard notation for the wave operator g°'^'Wa^ 

The first equation we need on our way towards Buchert averages is the Hamihonian con- 
straint or Goo component, obtained from contraction of the renormahzation - group improved 
Einstein equation (2.2) with u'^u", i.e. 



(2.6) 



where we denote by Kij the extrinsic - curvature tensor of the spacehke hypersurfaces that 
fohate the spacetime manifold and have taken T^j^^, of the form 14 1 



while 



]^Gl - GG,oo + ^(VpG)(V^G) - GUG 



(2.7) 



' A \ G J G G G G' 



(2.8) 



with A = g'^^'Vi'Vj the Laplacian operator (up to a sign). Hereafter, it is convenient to 
define 

ifjG = log G, (2.9) 

and bear in mind that the spatial part of the spacetime four -metric is, in our coordinates, 
the induced Riemannian three -metric hijdx^ (S> dxK The renormahzation - group improved 
00 component of the Einstein equation (2.2) reads therefore (we now write explicitly all 
traces, since the desired Buchert averages are for traces of the Einstein equations 14| ) 



(Trfs:)2 -TrfsT^ 

3 5 AG 

levrGp + 2A - -ij^ Q + e^Gfl + ^/i(grad?/^G, grad^/^c) - 



where, by definition. 



(grad^G)i = 



(2.10) 



(2.11) 



Since the left-hand side of (2.10) has the same functional form as in general relativity, the 
Buchert average (1.4) of (2.10) leads to (see appendix) 

2 



3(^) -M{Gp)d + \{^'^R)d-{^)d 

\0,D / ^ 

1 3 5 

-2<5d + {dipG,o)D - -^{^g,o)d + -(/i(gradV^G,grad?/^G))D 



'AG^ 
G , 



(2.12) 
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The second equation we need can be obtained by contracting (2.2) with the contravariant 
spacetime metric g^" . Of course, this contains also BPq already encoded, up to a sign, in 
(2.10), but the difference between the full scalar curvature and R\ is the new equation we 
need. We point out that, from (2.2) and (2.5), and exploiting the Arnowitt- Deser- Misner 
identity with unit lapse function and vanishing shift vector 



ot 



(2.13) 



one finds, by adding and subtracting TtK^ in (2.13), 



(^)i? = + (TtK? - TyK^ + 2TtK^ - 2—TtK 

L -I ot 



d 
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.G 



= 8nGp + 4A - -^^ + 39ijG,o + 3 ^ 

9 AG 
+ 2^(g^^'^^G',gradV'G) - 3-^, 



00 



(2.14) 



where we can exploit (2.10) and then take the Buchert average to obtain eventually (see 
appendix) 



\ dt /j 

2 2 d 



+ ^(/i(gradV'G,gradV'G))D - ^ (^^^ • 



(2.15) 



If the FLRW symmetry is irnposed, our equations (2.10) and (2.14) are in full agreement 
with (3.18a) and (3.18b) of 16|, respectively, i.e. {H being the Hubble parameter ^, and 
X = 1,0,-1 the curvature parameter, while bearing in mind that our pressure parameter 
vanishes) 



X 



1 

a? 3 



87r 

y 



A + ^Gp + ^G,o^ - 



a 



H^ + 2- + 



X 



a 



Otherwise they express corrections involving the spatial gradient of G. Note that 8tt{Gp)d 
can be eliminated in (2.15) with the help of (2.12), but this step is inessential. 

Last, but not least, we have to take the Buchert average of the on -shell consistency 
equation derived in [l6|, i.e. 

47r(G'r,W^G)B = (u^V^(GA))b, (2.16) 



where m'^ is the same vector used in (2.6). 



III. INTEGRABILITY CONDITION 



The integrabihty condition for our coupled system (2.12), (2.15) and (2.16) is obtained 
after re - expressing (2.12) and (2.15) in the form 



3 (^)' - 87r(Gp)^ + le^R)n - {A)d = -\Qd + (^i) 



an ^ 



where we have defined 



AG 



Fi = 9^/jg,o - |^G,o + ^^(grad?/^G, grad^/^c) ^ 



1 G 2 1 AG 

F-:,, = -OilJn n - n -\ iz- + -/i(grad'?/'G, gradi/jc) 



(3.1) 
(3.2) 

(3.3) 

(3.4) 



3 r^,u ru,u ■ G 3 3 G 

Following [18|], we now take the time derivative of (3.1) and then use again (3.1) and (3.2), 
i.e. (hereafter, all partial time derivatives acting on Buchert averages coincide with total 
time derivatives of such spatial averages) 

= 2— (s^ - 3 '^'^^ 



dt \ar)J 



an 



2— 

an 



4'k{Gp)d + {A)d + Qd + -{F2)d 



87t{Gp)d + \{^'^R)d - (A)d + - {Fi)d 



an 
Sir 
d 



-12n{Gp)o + + 1{^'^R)d + (^^2 - F,^ ^ 



{Gpe)n-3—{Gp)n 
an 



IOQd d 

D — -^T, r 77:\-f^l/D, 



2 dt dt 

where use has been made of the Buchert identity [l^ 



1^^^" -\dt/^ 



D 



(3.5) 



(3.6) 



with ip = Gp. Now both sides of (3.5) contain the term — 24:71 {hi)/ an) {Gp)^, leading 
eventually to the desired integrabihty condition 



at an ot 



D 



''D 



dt 



= IGtt 
+ 2 

where, from the definitions (3.3) and (3.4) 



'dt 

+ {Gpe)D 



+ 2|(A). 



+ — ((2Fi-3F2))b 
dt ud 



2Fi - 3F2 = 9ijG,o + li^G,o - 3^ + ^/i(gradV'G, grad^-c) 



AG 



(3.7) 



(3.8) 



IV. THE VARIABLE - G COSMIC QUINTET AND COSMIC EQUATIONS 

it is convenient to introduce the following dimensionless 



Inspired by the work in 
quantities: 



3if? 



D 



3i/? ' ^ 



D 



1^ 



(i^i) 



D 



'D '-'-'-'D 'J-'-'D 'J-'-'D "J-^D 

where the definition of has been suggested by the averaged Hamiltonian constraint 
(3.1), and hereafter = a£,/a£, is the effective Hubble parameter. In a homogeneous 
and isotropic universe with constant G and A, fi^ and reduce to the usual matter and 
cosmological constant density parameters, reduces to Q^, while Qq and Qq vanish, so 
that the standard FLRW cosmology is then recovered. Equations (4.1) represent therefore 
the definition of the cosmic quintet of density parameters for variable - G inhomogeneous 
models of the Brans- Dicke type. Such an interpretation is also confirmed on noting that, 
by virtue of (3.1), one finds 



(4.1) 



fi^ + + fig + fig + = 1 



G 



(4.2) 



Pursuing the analogy with FLRW cosmology, it is instructive to recast (3.1), (3.2) and (3.7) 
into an almost Friedmann form. To this aim, we look for an effective density and effective 
pressure such that (3.1) and (3.2) read as 



3i/2 = (A)^ + 87r(G)z)pffF, 



3^ = (A)^-47r(G)B(pfff + 3p, 



(4.3) 
(4.4) 
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By comparison with (3.1) and (3.2), we solve for and to find 

{G)nP^, = {Gp)o - -^e'R)D - ^Qd + (4.5) 
(G)nP?, = -^Qn + ^e'R)o - ^((Fi + 3F,))n. (4.6) 



V. ACCELERATING PATCHES AND STATIONARY MODELS 

From equation (3.2), inspired by the work in jl^, we see that the condition for an 
accelerating patch D of the universe reads as 

Qd > Ah{Gp)d - {A)d - 1{F2)d. (5.1) 
The definition of the density parameters in (4.1) can be used to re-express (5.1) in the form 

-n2> -nZ - -n^ - -^^^l (5.2) 



If the domain D is taken to be as large as our observable universe [18|, the Hamiltonian 
constraint (4.2), jointly with (5.2), yields 

|fif + ng>i-5ns-(||||^ + i)ng. (5.3) 

A. Stationary models 

In our approach, in order to evaluate the quantities occurring in the equations, stationar- 
ity of the whole universe model is assumed. Within this framework, the manifold S we refer 
to hereafter must be a finite-volume compact manifold, so that the global average makes 
sense. We find therefore the global stationarity conditions (cf jisl ) 

Qj: = AniGph - (A)s - ^{F^h, (5.4) 

= 127r (G'p)s - 6iJ| + 3(A)e + ((^^2 + 2Fi)^^ , (5.5) 

H^ = ^ = -, (5.6) 

as as 

where Qs and {^^^R)j: are now the global kinematical backreaction and the globally averaged 
three-dimensional spatial curvature, respectively, while the constant C can be obtained from 
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the Hamiltonian constraint evaluated at the initial time. Interestingly, if we now insert the 
global stationarity conditions (5.4)-(5.6) into the formulae (4.5) and (4.6) for effective density 
and effective pressure, we find some remarkable cancellations (i.e. vanishing coefficients of 
(G'p)e, (i^i)s, (i^2)s), leading to 



D 

eflf 



2' 



(5.7) 
(5.8) 



Hence we find the cosmic equation of state 

On taking the time derivative of Qs and (*^^-'-R)s in (5.4) and (5.5), and then using the 
identity (3.6) with ip = Gp,D = S, and again (5.4)-(5.6), we find 

d 



^ +3— )(Qs + (A)s + ^(F2)s) =47r 



dt 



(5.10) 



ot as at 



{GpO), 



(5.11) 



The complete integral of (5.10) is given by the complete integral of the homogeneous 



equation 



m 



plus a particular integral of the full equation, so that we can write (the initial value of the 
left-hand side below being denoted by Js(ij)) 



+ 47r f\{t,t') 



' d_ 
dt 



-{Gp)^ Jt') + {Gpeut') 



dt' 



(5.12) 



where 'y{t,t') is the Green function of the operator |^ + 3^^^^. Since we are here assuming 
as = C, we are actually dealing with the first-order operator ^ + j^^^^, with aj = aj]{ti). 
The desired Green function solves the equation 



d 



+ 



3C 



dt (Ct + ai)^ 



i{t,t') = oyt^t', 



(5.13) 
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and suffers a jump aX t = t' given by 19 1 



lim 7(t,t') - lim 7(t,t') = 1. 
t^t'+ t^t'- 

It therefore reads as (hereafter G is the step function) 

7(t,t') = {Ct + ai)-^^\Q{t ~ t')h^{t') + Q{t' - t)h2{t') 



where 



i.e. one finds, Vt 7^ t', 



3C 



(5.14) 



(5.15) 



(5.16) 



(5.17) 



VI. A SOLUTION FORMULA FOR THE SPATIALLY AVERAGED SCALAR 
CURVATURE 



Note now that, by virtue of (5.10), the equation (5.11) can be recast in the form 
- ( - 3Q^ - 6(A)s - im + 2i^i))s) = -G^M = 12^ 



dt 



dt 



3 ' 



(6.1) 



where a remarkable cancellation of coefficients for ^Qt, has occurred. On denoting by x ^'^ 
integration constant, and bearing in mind that cie(^) = Ct + ai, equation (6.1) is solved by 

= 3Qj: + 6(A)s + im + + X + [ ^^^^^dt\ (6.2) 

where is given by the solution formula (5.12). The coupling between spatially averaged 
scalar curvature and backreaction is therefore found to survive. 



VII. A PARTICULAR ASSUMPTION ON THE AVERAGED EQUATIONS 

In order to work out a general solution of the cosmic equations, we should know how 
matter is spatially distributed and assume a metric for the inhomogenous universe to com- 
pute the averaged quantities. Moreover, an ansatz should be made for the spatial variation 
of G, possibly consistent with the renormalization-group equations, so that Fi and F2 can 
be evaluated. Since such an information is lacking by virtue of severe technical difficul- 
ties in deriving G{x,t), we can only look for a particular solution under some reasonable 
assumptions. 
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A first step along this road can be done by differentiating both sides of Eq. (4.3) with 
respect to the cosmic time t, thus yelding 

= ^^^^ + ^"^^^^ + '^^^ ) + ^ [~dr + ^"^^^^r j • 

In a FLRW universe, the effective fluid should be replaced by the standard source (matter 
and radiation) term. In this case, the standard continuity equation holds and the above 
relation identically reduces to the Raychaudhuri equation. Let us assume that our effective 
fluid can still satisfy, to a first approximation, the standard continuity equation, i.e. 

^ + 3//z,(pfff+pfff)=0. (7.1) 

In classical general relativity, this equation holds by construction of the effective equations. 
In our approach, consistency of the theory seems to demand that Eq. (7.1) should hold, 
but for example we cannot yet say what would one expect if (7.1) had source terms. To 
study (tiny) departures from (7.1), one might try to assume analyticity of G(x, t), A(x, t) 
and hence set up a perturbative scheme for the evaluation of all averaged equations. This 
is an important topic that deserves attention in a separate paper. 

On solving with respect to {l/HD)dp^fi/dt, the above relation becomes 

By equating to Eq. (4.4), we get the remarkable relation 

^^ + 8vrpe,^^ = 0, (7.2) 
which, for a homogenous and isotropic universe, reduces to 

A + SvrpG = 0, 



which has already been obtained in the literature |ll|. 

It is worth wondering whether Eq. (17. ip can lead to a constraint also on the functions Fi 
and F2. For this purpose, let us first note that, by virtue of the definitions of the effective 
fluid energy density and pressure, one has 

d{G) 

^ '^~dr~~dt 16^^r~16^ di ^~dt 
(G')D(pfff = {Gp)d - I^Qd - i^{'R)d + - ^((Fi + m)D. 
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Multiplying by {G)d and inserting these relations, the continuity equation (17.11) for the 
effective fluid becomes 



D 



dt 



+ QHdQd 



dt 



D 



D 



levr 



( d{Gp) 



D 



+ 2 



V dt 



+ 3Hd{Gp)d 



dt 



dt ■ 

Comparing this relation with the integrability condition (3.7) and making use of the Buchert 
identity (3.6) with ip = Gp and {9)£) = SH^, we finally get 



d{A) 



D 



+ 8ivp^s^{G)d = Hd(HF^)d - 2((Fi + ?,F2))d - ((2Fi - 3^2)),^ 



dt dt 

so that, since the left-hand side vanishes because of Eq. (17.21) . we obtain 

3(Fi)^ - 2((Fi + 3^2))^ - ((2Fi - 3F2))b = -((Fi + 3^2))^ = 0, 



(7.3) 



having made use of the linearity of the mean operator. Summarizing, we can therefore 
conclude that the conservation of the effective fluid implies two remarkable relations. The 
first one, Eq. (17.21) is the counterpart for an inhomogenous model of the renormalization- 
group condition on the G and A time-derivatives. On the contrary, Eq. (17.31) is an interesting 
new constraint closely relating the spatial variation of the two auxiliary functions Fi and 
F2. Indeed, whatever is the exact spatial variation of G, the condition (-Fi)d = ~^{F2)d 
must hold. It is worth stressing, however, that this constraint only refers to the averaged 
quantities, so that it is impossible to infer any information on the spatial variation of G. 

In order to gain further insight, it is interesting to consider the case of a nearly homoge- 
nous and isotropic universe. As a first step, we introduce three auxiliary effective fluids with 
energy densities 



n 



{Gp)d 3 {F2)d 



D 



{G)d ' 



1 Q 



D 



IQtt (G) 



D 



(7.4) 



D_ 1 {'•'^R)d , 1 {{Fi + 3F2))d 
Pr ~- 



IQtt (G) 



D 



ill 



(G) 



D 
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and pressure given as 




0, 



< 




1 Qd 
167r {G)d 



(7.5) 




1 ((3)/?) 



D 



1 ((Fi + 3F2))^ 
247r (G)!, 



487r (G)^ 



so that the equations of state are {w^j, Wq, w^) = (0, 1, —1/3). As can be straightforwardly 
checked, one has 



which shows that the effective fluid may be viewed as consisting of three distinct components. 
The first one, p^, has zero pressure and reduces to the standard matter term when we turn 
off the spatial variation of G (so that Fi = F2 = 0) and replace averaged quantities with 
standard ones in a FLRW universe. We can therefore think of it as the matter term of 
our model although, since we do not know anything about Fi and F2, we cannot rule out 
a priori that {—3/8'k){{F2)d/{G)d) overcomes {Gp)^/ {G)d^ thus giving rise to a negative 
energy density. The second component, with energy density p^, has the same equation of 
state as the stiff matter. However, should the backreaction term be positive, its energy 
density is negative thus leading to a negative pressure acting as a variable cosmological 
constant (eventually driving cosmic speed up). Finally, the term with energy density p^ 
has a negative equation of state, but we cannot say whether it acts as a speeding up factor. 
Indeed, we do not know whether the curvature term overcomes the other term giving a 
positive or negative energy density. Note that, should we assume the conservation of the 
effective fluid, then Eq. (17.31) follows and p^ is positive-deflnite. However, to obtain results 
of general nature, we prefer not to assume here a priori the validity of Eq. ( 17. ip . In other 
words, at this stage, we are considering effective fluids which make it unnecessary to consider 



On denoting with w^^ = Pes/PeS the effective fluid equation of state, and inverting Eqs. 



PcS = Pm + Pq+Pr^ 



< 



(7.6) 



[PeS=PM+PQ+PR: 



Eq. (7.1). 
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(rrejl . we get 



^pfl 



-l + (l + 3u;, 



D \PeS 



(7.7) 



D \Pea 
''^P^M 



which are fully general. Let us now assume that the universe is nearly homogenous and 
isotropic as it is expected to be, for instance, in its infancy when perturbations had still to 
grow. In this case, we can write 

Gp - (3/87r)F2 = Gjvpr/''''^[1 + ^Mit,X% (7.8) 

G = GN[l + AG{t,X% (7.9) 

where Gat is the Newtonian gravitational constant, p^/"^*^ the dust matter energy density in 
a FLRW universe, and AAf(t, X*) and AG(t, X*) two unknown functions accounting for the 
small perturbations induced by the inhomogeneities. It is worth noting that Eqs. (17. Sp and 
(17. 9p may still be formally written even if deviations from the FLRW universe are severe, 
the only difference being that Ajvf(t, X*) and AdtjX^) are no longer much smaller than 1. 
On averaging Eqs. (17.81) and (17. 9p . the energy density of the matter -like term reads as 

(7.10) 



^ _ {Gp-{3/8n)F2)D 
Pm — 



D 



D 

Pm ^ ■ 



1 + (Ag)d 



while Eqs. (4.5) and the average of (17. 9p give 



3Hl 



(A) 



D 



8716-^(1 + (Ag)d)' 
On inserting Eqs. (17.101) and (17.111) into Eqs. (17. 7p . we get 



(7.11) 



^ ^ 1 / 1 + {Am)d 



J, _3 fl + {Am)d 
pR - 4 



-1 



1 + (Ag) 



D 



1 + {Am)d (87rG^pi;/^^) 



1 + {Ag)d {SHl - {A)n) 



(l + 3w, 



eS) 



'eff) 



^FLRW 
Pm ) 



FLRW 

Pm 



(7.12) 



1 + (Am)i. (87rG^pf/^^) ^ 
We can now work out an interesting property of these two effective fluids by assuming a 
phenomenological ansatz for the term 3H^ — (A)d and for the effective fluid equation of 
state w^. It is reasonable to assume that, for very large 



FLRW / 
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Equation (17.121) shows that, in this hmit, both pq and vanish so that the dust matter 
FLRW case is recovered in the early universe whatever is the exact shape of the correction 
functions Am and Aq- In other words, we find that the early universe tends to a FLRW 
model, but keeps track of the original inhomogeneities through the two effective fluids with 
density pg and p^. Actually, depending on the functional expression adopted for Hd, it is 
also possible that (A^) vanishes identically at all redshifts. We stress that the constraint 
{Am)d{z) = does not imply that the universe is homogenous, but only that the matter 
inhomogeneities average out to zero. As is clear from Eqs. (17.121) . in such a case, the two 
fictitious fluids pg and p|j still contribute to the cosmic dynamics. In particular, should 
Hd lead to cosmic acceleration, we can therefore argue that these two fluids arising from 
inhomogeneities drive the cosmic speed up in a matter dominated universe. 



VIII. CONCLUDING REMARKS AND OPEN PROBLEMS 



Our original equations (2.12), (2.15), (2.16), (3.7), (4.5), (4.6), (5.2), (5.4), (5.5), (5.10), 



(5.12), (6^ 



averages 



14 



provide a simple but nontrivial application of the Buchert method for spatial 



18| to renormalizat ion-group improved action functionals with variable G and A 
of the Brans-Dicke type [16] . It now remains to be seen whether our averaged equations agree 



with the qualitative picture in general relativity 18|, l20l], according to which backreaction 
effects point to a global instability of the standard cosmological model, with exact solutions 
and perturbative results modeling this instability lying in the right sector to account for 



dark energy from inhomogeneities 



21 



22j (for a critical view, see however the work in |23| ) . 



The definition of the effective fluids in section 7 makes it clear what is the role played 
by the different inhomogeneity terms introduced by the averaging procedure. In the early 
universe, only the p^j term survives and reduces to the standard dust matter, so that the 
usual FLRW case is recovered. It is however intriguing to note that, in the late epochs, the 
effective matter fluid energy density p^ is increased with respect to the standard one because 
of the {—S/87!'){F2)d/{G)d term. As such, one could speculate that this latter component 
mimics an effective dark matter component, while {Gp)d/{G)d accounts for the baryons. 
On the other hand, the two additional fluids pg and p^ can both provide a negative pressure 
so that they can drive accelerated expansion. The resulting qualitative picture is that of a 
universe consisting of baryons only, while inhomogeneities average out to give rise to the full 



17 



dark-side phenomenology. Needless to say, such a qualitative picture must be substantiated 
by a detailed comparison with the data on both the background expansion and the growth of 
structure, and hence we cannot draw any definite conclusion at the moment. An alternative 
possibility, from this point of view, is to compare our effective fiuids' picture with the recent 
literature on the morphon field 



24 



25( 1 in order to work out fruitful analogies. 



We may be criticized for having used the renormalization-group approach only as a moti- 
vation, without ever exploiting it for the explicit evaluation of the Buchert averages involving 
G and A. On the other hand, as far as we know, no-one has succeeded so far in obtaining 
G{x,t) and A{x,t) from the exact integration of the renormalization-group equation, and 
previous attention had always focused on FLRW models where, by symmetry, G and A can 
only depend on the time variable. Without an explicit knowledge of the desired G{x,t) and 
A{x,t), we cannot yet provide examples of our averaging procedure. However, as we have 
stressed after Eq. (7.1), a promising way out might be obtained by studying a power-series 
expansion of G and A, which is independent of any fixed-point assumption. We hope to be 
able to return to this point in a separate paper. 

Another interesting issue is whether our averaged equations with variable Newton pa- 



rameter can be relevant for the theoretical scheme proposed in 



26 



21, 



28j as yet another 



alternative to dark energy. It would be also quite important to repeat our analysis with the 



help of the covariant technique developed in 



29|. 
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APPENDIX A: BASIC IDENTITIES FOR THE BUCHERT AVERAGES 



Relying upon 



14( 1 . we express the extrinsic-curvature tensor Kij of the spacelike hyper- 



surfaces foliating the spacetime manifold in the form 

e 



— cr, 



(Al) 
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where the shear tensor is symmetric and traceless. Hence one finds the simple but funda- 
mental identities (hereafter = ^aija^^) 



3 ' 

-((Trir)2-Trir2) = ^_^2 



(A2) 
(A3) 
(A4) 



2^' ' '3 
When we evaluate the Buchert average of the Hamiltonian constraint (2.10), we exploit the 
identity (1.7) and the definition (1.8) to write 



1 

3* 



Qd l 
2 3 



+ ^{0)1 



Qd 



+ 3 



an 



(A5) 



Along the same lines, the first equality in (2.15) is obtained by virtue of 



'ttK' - ^Trir^^ = 2{a')o + ^ [{9')n - {9)l\ + 1(9)1 + 



'd9\ 



where 11 



'd9\ _ d 



{9)d-{9')d+{9)1. 



(A6) 



(A7) 



The identities (1.7), (A6) and (A7) lead eventually to the first equality in (2.15). 
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